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Abstract 

The yields of (anti-)protons were measured by the NA49 Collaboration in centrality selected 
Pb+Pb collisions at 40^4 GeV and 158^4 GeV. Particle identification was obtained in the laboratory 
momentum range from 5 to 63 GeV/c by the measurement of the energy loss dE/dx in the TPC 
detector gas. The corresponding rapidity coverage extends 1.6 units from mid-rapidity into the 
forward hemisphere. Transverse mass spectra, the rapidity dependences of the average transverse 
mass, and rapidity density distributions were studied as a function of collision centrality. The values 
of the average transverse mass as well as the midrapidity yields of protons when normalized to the 
number of wounded nucleons show only modest centrality dependences. In contrast, the shape of 
the rapidity distribution changes significantly with collision centrality, especially at 40A GeV. The 
experimental results are compared to calculations of the HSD and UrQMD transport models. 
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I. INTRODUCTION 



It is generally accepted that heavy ion collisions at ultra-relativistic energies result in 
a fireball of matter with high density and temperature. Such conditions prevail, when 
the incoming nucleons deposit enough of their kinetic energy in the reaction zone. Little is 
known about this stopping process and about possible differences between the stopping of the 
incident nucleons in elementary nucleon-nucleon interactions and nucleus-nucleus collisions. 
In a simple microscopic picture each nucleon interacts only once in elementary interactions, 
whereas in central nuclear collisions multiple interactions prevail. There are three methods 
to study nuclear stopping. In proton-nucleus (p+A) collisions the difference between the 
c.m.-energies of the incident and the most forward going nucleon is a good measure of the 
stopped energy. This type of analysis has been pioneered in reference [1]. Another approach 
to study effects due to multiple projectile nucleon collisions employs symmetric collisions 
of nuclei of different size. The resulting distributions of participating nucleons in terms of 
longitudinal and transverse momenta may vary with the size of the incident nuclei. Such 
differences reflect the change of stopping power as a function of system size. Finally, the 
effective system size can be changed by varying the impact parameter in collisions between 
heavy nuclei. Here again the momentum distribution of participating nucleons may be 
different in central and peripheral collisions as a consequence of changes in the stopping 
behavior. 

Energy loss in central collisions was studied as a function of beam energy at AGS in 
Brookhaven [2], at CERN SPS in Geneva [3] and at RHIC in Brookhaven [4]. It was 
found that the rapidity distribution changes from a convex to a concave form from AGS to 
RHIC via SPS energies. At AGS energies the energy loss was also studied as a function of 
collision centrality [2]. A strong centrality dependence was observed. This topic has also 
been addressed by NA49 in two previous publications. The first one presents the net-proton 
distribution for central Pb+Pb collisions at 158 A GeV [3]. The second one describes the 
measurement of midrapidity proton and antiproton yields in Pb+Pb collisions at various 
energies from 20A to 158 A GeV as well as for different centralities at 158A GeV [5]. At 
RHIC the PHOBOS collaboration has published the centrality dependence of the rapidity 
density at midrapidity and the transverse mass spectra of net protons [6] at y/s^ = 62.4 
GeV and charged particle pseudorapidity distributions [7, 8] at y/s NN = 19.6 GeV, 62.4 GeV, 



and 130 GeV. Data on centrality dependent particle production at y = and y = 1 in 
Au+Au collisions at ^ l /s^ = 200 GeV [9] are available from the BRAHMS collaboration. 
The STAR collaboration has published Au+Au data on the centrality dependence of proton 
and antiproton production for \y\ < 0.5 at y/s NN = 130 GeV [10]. 

In this paper we present the centrality dependence of proton and antiproton transverse 
mass and rapidity distributions in Pb+Pb collisions at 40A GeV and 158^4 GeV as obtained 
with the NA49 detector [11]. The phase space coverage extends in center-of-mass rapidity y 
from mid-rapidity 1.6 units into the forward hemisphere and ranges from zero to 2.0 GeV/ c in 
transverse momentum pp- Net-proton distributions at 158 A GeV are obtained from those 
of all protons by subtracting the distributions of antiprotons. This analysis supplements the 
NA49 data on n, K, <p, A production as function of beam energy at the CERN SPS [12-14]. 
Pion and kaon production in Pb+Pb collisions at A0A GeV and 158A GeV is analyzed as 
function of centrality with the same method and is the subject of a separate paper [15]. 

II. EXPERIMENTAL SETUP AND DATA SETS 

The NA49 detector is a large acceptance hadron spectrometer at the CERN SPS [16]. The 
main components are four large time projection chambers (TPCs) and two super-conducting 
dipole magnets with aim vertical gap, aligned in a row, and a total bending power of 9 Tm. 
Two 2 m long TPCs (VTPCs) inside the magnets each with 72 pad-rows along the beam 
direction allow for precise tracking, momentum determination, vertex reconstruction, and 
particle identification (PID) by the measurement of the energy loss (dE/dx) in the detector 
gas. The other two TPCs (MTPCs)) have large dimensions (4m x 4m x 1.2m, 90 pad-rows) 
and provide additional momentum resolution for high momentum particles as well as PID by 
dE/dx measurement with a resolution of around 4%. Two time-of-flight scintillator arrays 
of 891 pixels each, situated just behind the MTPCs symmetrically on either side of the beam 
axis, supplement particle identification in the momentum range from 1 to 10 GeV/c. A Veto 
Calorimeter (VCAL), which is placed further downstream along the beam and covers the 
projectile spectator phase space region, is used to select event centrality. The NA49 detector 
is described in detail in reference [16]. The Pb beam had a typical intensity of 10 4 ions/s 
and impinged on a target Pb foil with a (areal) density of 224 mg/cm 2 . It passed through 
a quartz Cherenkov detector from which the start signal for the time-of-flight measurement 



TABLE I: Cross section fractions in % centrality, average numbers of wounded nucleons (iV w ), 
and numbers of analyzed events for the five centrality classes at 40A GeV and 158^4 GeV. Only 
systematic errors are quoted. 



centrality class 


centrality [ % ] 


(N w ) 


analyzed events 


AOA GeV 


CO 


0-5 


351 ± 3 


13034 


CI 


5-12.5 


290 ± 4 


22971 


C2 


12.5-23.5 


210 ± 6 


34035 


C3 


23.5-33.5 


142 ± 8 


32668 


C4 


33.5-43.5 


93 ± 7 


32071 


158 A GeV 


CO 


0-5 


352 ± 3 


15306 


CI 


5-12.5 


281 ± 4 


23548 


C2 


12.5-23.5 


196 ± 6 


37053 


C3 


23.5-33.5 


128 ± 8 


34554 


C4 


33.5-43.5 


85 ± 7 


34583 



was obtained, and three stations of multi-wire proportional chambers which measured the 
trajectories of individual beam particles. A minimum bias trigger was derived from the signal 
of a gas Cherenkov device right behind the target. Only interactions which reduce the beam 
charge and thus the signal seen by this detector by at least 10% are accepted. The interaction 
cross section thus defined is 5.7 b. The contamination by background events remaining 
after cuts on vertex position and quality amounts to less than 5% for the most peripheral 
collisions and is negligible for near-central collisions (see [17]). The resulting event ensemble 
was divided into five centrality classes CO, CI, C2, C3, and C4 (see Table I and [17]). The 
centrality selection is based on the forward going energy of projectile spectators as measured 
in VCAL. The number of projectile spectator nucleons and the average number of interacting 
(wounded) nucleons (N w ) were calculated from the selected cross section fractions using the 
Glauber approach. The track finding efficiency and dE/dx resolution were optimized by 
track quality criteria. To be accepted, a track must have at least 50 (out of a maximum of 



90) potential points in the MTPCs and have at least 5 measured and 10 potential points 
in one of the VTPCs. Finally, tracks were required to have azimuthal angles within +/- 30 
degrees with respect to the bending plane in order to minimize reconstruction inefficiencies 
and to optimize the accuracy of the dE/dx measurements. 



III. ANALYSIS METHOD 
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FIG. 1: a) The most probable specific energy loss e for particle species 7r + , K + und p as a function 
of total momentum p. The points show results extracted from data whereas the lines indicate the 
Bethe-Bloch parametrizations optimized for the NA49 measurement, b) The ratios of the Bethe- 
Bloch parametrizations (ej) to from data are shown as function of laboratory momentum. The 
index i stands for p and p. 

The protons are identified by the measurement of their specific energy loss in the rela- 
tivistic rise region. As an appropriate measure of dE/dx we calculated for each track the 
truncated mean of the distribution of charges measured in each pad row of the MTPCs. 
Their raw yields were extracted by fitting the function F to the dE/dx distribution of all 
positively charged particles in narrow bins of total momentum p and transverse momentum 
Pt [18]. Antiproton yields are determined from the dE/dx distributions of all negatively 
charged particles with the same method at 158A GeV. At 40A GeV the antiproton statistics 
was very low and did not allow reliable extraction of yields. The shape of F is assumed to 
be the sum of Gaussians. Their parameters depend on the particle masses and the mea- 
sured track lengths. We modified the Gaussian functions by means of an extra asymmetry 



parameter to account for tails of the Landau distributions which are still present even after 
truncation. 

The function F reads 



\ dx J i=dJ?K,«,e Si U l I V27T^,Z 



exp 



1 ( ej(p) -etipY 



(1) 



Here ^ (abreviated e) is the measure of specific ionization, and e is its most probable 
value. The other parameters of the function are: 

• Af the raw yield of the particle % under consideration in a given phase space bin. 

• nf. the number of tracks in a given track length interval /. The second sum together 
with the normalization J2i n i forms the weighted average of the track ensembles in 
each phase space interval. 

• the width for the asymmetric Gaussian of particle type i in length interval /. 

• 5: the asymmetry parameter. 
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FIG. 2: Raw yields of protons at 40A GeV (a), protons (b) and antiprotons (c) at 158^4 GeV are 
shown as function of center-of-mass rapidity y and for centrality class C2. 



The amplitude parameters Aj(p,px) were determined by a maximum likelihood fit to 
the dE/dx histogram in each p, p^ bin. The peak positions e"j(p) are considered to be 
Pt independent and were determined by fits to the integrated distributions in all p 
bins. The fitted peak positions Si(p) for ir + , K + and p are compared to a parametrization 
of the momentum dependence of the most probable energy loss determined for the NA49 
detector in Fig. la. The relative differences between fitted (anti-)proton positions and the 



Bethe-Bloch parameterization are smaller than 5% (see Fig. lb) . We have verified that 
the peak positions do not depend on centrality. The widths of the Gaussians (cr^) depend 
on the particle type i and on the track length / according to cr^ = a ■ (ei/e 7r )°(l/\/T). 
The exponent a (= 0.625) was extracted from simultaneous fits to m? distributions from 
TOF and to dE/dx distributions from the TPCs [18]. The momentum averaged widths (a) 
were determined for each centrality bin by fits to the data in the whole p range. These 
momentum averaged widths (a) turned out to be approximately 4 % for each centrality bin. 
5 was studied by fits with the asymmetric Gaussians. Since it did not show any significant 
variation with centrality, total and transverse momentum [11], this parameter was fixed to a 
constant value (0.071). The total number of fit parameters for each p, p^ bin is 9. The raw 
yields were transformed from a fine grid in log(p), pr to a coarser grid in y, px- Examples 
of the transformed results for centrality class C2 are shown in Fig. 2. 

A. Acceptance and inefficiency 

The raw particle yields have to be corrected for losses due to tracks which do not pass 
through the detectors or which do not fulfill the acceptance criteria (acceptance losses) 
and tracks which are not properly reconstructed (efficiency losses). The acceptance was 
calculated by generating a sample of (anti-)protons with flat distributions in transverse mo- 
mentum and rapidity. The generated particles were propagated through the detectors (and 
the magnetic field) using the programs provided in the GEANT 3.21 [19] package. Along 
the resulting trajectories realistic detector signals were generated and processed in exactly 
the same way as the experimental data. The ratio of generated to accepted particle tracks 
in each y, pr bin is the acceptance correction factor. In addition to the well-defined correc- 
tions necessary to correct for the limited and constrained acceptance, the raw spectra may 
be subject to losses due to detector occupancy and thus centrality dependent inefficiencies. 
These losses were minimized by restricting the analysis to tracks with azimuthal emission 
angles within ± 30 degrees with respect to the bending plane and in the bending direction. 
The remaining losses were determined by the following procedure: Ten GEANT generated 
(anti-)proton tracks and their signals were embedded into raw data of real events. Only 
those tracks are embedded which pass all acceptance criteria. These modified events are 
reconstructed with the standard reconstruction chain. The ratio of all generated tracks to 



those reconstructed constitutes the applied correction factor for reconstruction inefficiencies. 
The resulting efficiencies vary with centrality by less than 5% and are over 95% in most y, pr 
bins. In the further analysis the bin size in rapidity of the antiproton spectra (see Fig. 2c) 
was doubled in order to reduce the statistical errors. The one in was also doubled at 
AOA GeV. 

B. Feed-down corrections 

The measured (anti-)proton yield contains (anti-)protons from weak hyperon decays, 
namely of the A (A), the S° (E°), and the S + The (anti-)proton contribution from 

S + decays is determined by scaling the estimated feed-down correction coming from 

A (A) by the ratio A / S + (A derived from a statistical hadron gas model [20] thus 

assuming the same phase space distributions for A (A) and S + hyperons. To determine 
the feed-down correction from A (A) in each measured phase space bin, we used a similar 
procedure as for the efficiency calculation. A (A) were generated according to distributions 
measured by NA49 [14, 21] and embedded into real events. Those (anti-)protons from 
embedded A (A) decays, which are reconstructed and accepted by the track selection cuts as 
such, have to be subtracted from the (anti-)proton yield after proper normalization per event 
(see Eq. (2)). In reference [14, 21] the measured A (A) include the A (A) from electromagnetic 
decays of S° und S° . The feed-down correction is given by: 



embedded A (A) decays, N^^Jy , p?) is the number of simulated A (A) and Y tot {y,pi) the 
multiplicity of A (A) in the phase space interval from reference [14, 21] scaled for the contri- 
butions from S + (S~). The lowest and the two highest rapidity bins in the 40 A GeV data 
(see Fig. 2a) were removed due to low statistics in too many p^ bins. 

C. Systematic Errors 




Y to t{y,PT) 



(2) 



where N^°^ 



(|/,Pt) is the average number of reconstructed (anti-)protons per event from 



One of the sources for systematic errors is the uncertainty in the procedure of unfolding 
the dE/dx distributions. Studies of the sensitivity to small changes of the fit parameters 



showed that the systematic errors are of the same order as the statistical errors. The largest 
contribution to the systematic error was traced to the deviations of the peak positions from 
the predicted values (see Fig. 1). They reach up to 4% in the yields in the low momentum 
bins. The next source of systematic errors is the uncertainty related to acceptance, ineffi- 
ciencies and feed-down correction calculations. To estimate the magnitude of these errors 
cut parameters were varied such that significantly different correction factors were obtained. 
The resulting final results varied only within the statistical errors. Overall it appears that 
the systematic uncertainties resulting from the correction procedures are of the order of 3% 
or less. 

The determination of the rapidity density distributions required extrapolations of the 
transverse momentum distributions. Their contribution is negligible for most of the rapidity 
bins except the rapidity bins y > 1.2, where the extrapolation factors reached values of 1.3. 
Single and double exponentials were used to describe the shape of the transverse momentum 
distributions in this region. The resulting differences of pr-integrated yields are smaller than 
2%. Finally, the results from the dE/dx analysis in this paper and the TOF-results [5] agree 
within 5 % in the common acceptance region at all centralities. We conclude from these 
studies that each data point carries a systematic uncertainty of approximately 7 %. 

IV. RESULTS AND COMPARISON WITH MODELS 
A. Transverse momentum spectra 

Transverse momentum distributions of (anti-)protons were determined in 10 (11) bins 
of rapidity at 40 A GeV (158^4 GeV) as shown in Fig. 2. The range covered in trans- 
verse momentum extends from p^ equal zero to 0.8 GeV/c at high rapidity and up to 
1.5 GeV/c (2 GeV/c) for 0.0 < y < 1.0 at A0A GeV (158 A GeV). Whenever necessary 
the spectra were extrapolated to 2 GeV/c (neglecting contributions at higher trans- 
verse momenta which have been accounted for in the systematic uncertainties) by taking 
the mean of the single and double exponential fits to the measured invariant spectra. At 
A0A GeV the extrapolation of the two bins around midrapidity were performed with the 
functional form obtained from the adjacent rapidity bin because of low statistics in the 
Pr distributions. Using theses extrapolations px integrated yields (dn/dy) and mean val- 
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FIG. 3: Invariant mp spectra of protons at 40^4 GeV (a) and 158A GeV (b) as well as antiprotons 
at 158^4 GeV (c) for five centrality bins in Pb+Pb collisions. The respective center-of-mass rapidity 
intervals are - 0.02 < y < 0.18, - 0.12 < y < 0.08, - 0.12 < y < 0.08. The data 
at different centralities are scaled down by the factors indicated in the figures. The new NA49 
measurements (full symbols) are compared, whenever available, to results of an earlier analysis 
using TOF information for particle identification [5] (open symbols). Only statistical errors are 
shown if larger than the symbol size. For the systematic errors see section IIIC. 

ues ((m-p) — mo) as function of rapidity were calculated (see below). Midrapidity invariant 
mx spectra (l/mr d 2 n/(dmx dy)) at A0A GeV and 158A GeV are compared in Fig. 3 with 
the results of a TOF based analysis published previously [5]. We find agreement within 
errors except for some of the two lowest tut points with deviations up to 20% in the 
158A GeV data. These differences can be traced to the feed-down corrections, which relied 
on different (older) parameterizations used for the MC input of (anti-)lambda phase space 
distributions in reference [5], but are based on recent measured data in this analysis (see 
above). However, this uncertainty has only a small effect on the value of dn/dy. 

Since the shapes of all invariant spectra deviate significantly from single exponentials, 
we choose (m-p) — rno instead of the inverse slope parameter of the transverse mass spectra 
to study the transverse activity as function of rapidity. Fig. 4 shows (m-p) — m of protons 
at 40 A GeV (upper row) and 158 A GeV (middle row), and antiprotons at 158 A GeV (lower 
row) as function of cms rapidity. In the latter the bin size in rapidity was doubled to reduce 
the statistical errors on the data points. The (mp) — ttlq values near midrapidity are plotted 
in Fig. 5 as function of centrality. Also shown are the data points from the TOF analysis [5] 
and recent results from p+p interactions at the same energy [22]. A clear increase of 



(mx) — frio by roughly 30% is observed when comparing the values obtained in the most 
peripheral with the most central event sample. The increase is close to a factor of two when 
the results on p+p interactions are taken as reference. Similar observations have been made 
by NA49 for data on hyperons [21]. 



B. Rapidity spectra 

The rapidity densities dn/dy in each rapidity interval were obtained by summing the 
yields in the measured px interval and the integrals of the extrapolation function above 
the highest measured px bin. The additive correction was calculated as the average of 
the single and the double exponentials in the unmeasured px region. The contributions of 
the extrapolation to dn/dy are mostly of order of 3 % (6 %) for protons and below 5 % for 
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FIG. 4: (mx) — mo for (anti-)protons as function of rapidity for five centrality intervals in Pb+Pb 
collisions. We show in the upper row (a) protons at 40^4 GeV , in the middle row (b) protons at 
158A GeV, and in the lower row (c) antiprotons at 158^4 GeV. The open symbols are obtained by 
reflection at mid-rapidity. Asterisks refer to results of an earlier analysis using TOF information 
for particle identification [5]. Only statistical errors are shown (if larger than the symbol size). For 
the systematic errors see section IIIC. 
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FIG. 5: (mx) — mo near mid-rapidity as function of (N w ) in Pb+Pb collisions. The leftmost panel 
(a) shows the proton data at 40^4 GeV, the middle panel (b) the proton data at 158A GeV, and the 
right panel (c) the antiproton data at 158^4 GeV. The results of this analysis (labelled dE/dx) are 
compared to previously published Pb+Pb data (labelled TOF) [5]. The ordinate has a suppressed 
zero. Only statistical errors are shown (if larger than the symbol size). For the systematic errors 
see section IIIC. 

antiprotons at 158A GeV (40A GeV) and contribute negligibly to the errors of the integrals. 
These extrapolations reach up to 25 % close to the limits of the accepted rapidity region. 

The rapidity spectra of protons at A0A GeV (left) and net protons at 158A GeV (right) 
are presented in Fig. 6 for five different centrality selections (CO - C4). We included the 
dn/dy values obtained from the TOF analysis (asterisks) published earlier by NA49 [5] 
for comparison. At both energies the yields increase with centrality. No change in shape is 
apparent at 158A GeV beam energy, whereas at 40A GeV the form of the dn/dy distribution 
evolves from a parabolic shape near to midrapidity (\y\ < 1) in semi-peripheral to a double 
hump structure in central collisions. The dn/dy values of protons at 40 A GeV are given in 
Table II, those for protons at 158A GeV in Table III and for antiprotons in Table IV. 

The trends in the evolution of the rapidity distributions are seen best, when the spectra 
are divided by the number of wounded nucleons as shown in Fig. 7. The proton spectra at 
158A GeV (middle panel) change little. The normalized yields (at midrapidity) increase by 
roughly 15% from 0.065 to 0.075 when going from semi-peripheral to central collisions. At 
AO A GeV the scaled yields at midrapidity increase by 25% from 0.085 to 0.11 from semi- 
peripheral to central collisions. As to the shape it seems that with decreasing centrality 
additional protons populate the region \y\ < (y max — 1-2). The shape of the antiproton 
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FIG. 6: dn/dy of (net) protons for five different centralities. We show in the left panel (a) the 
proton distributions at 40^4 GeV and in the right panel (b) the difference between proton and 
antiproton spectra (i.e. the net proton spectra) at 158A GeV. The open symbols are obtained by 
reflection at mid-rapidity. Also shown are the results from an earlier analysis labelled TOF [5]. 
Only statistical errors are shown (if larger than the symbol size). For the systematic errors see 
section IIIC. 

distributions at 158A GeV resembles a Gaussian and does not change with centrality. The 
integrated and mid-rapidity yields normalized by (N w ) decrease by 20% with increasing 
centrality. 

We present the dependences on centrality of the normalized antiproton multiplicity and 
of the width of their rapidity distribution in Fig. 8 together with NA49 data from elementary 
p+p interactions [22] both at 158A GeV. The mean multiplicities (total yield per event) 
of antiprotons for different centralities were calculated by integrating the measured rapidity 
spectra and by extrapolating into the unmeasured regions assuming a Gaussian shape (see 
Fig. 7). The magnitude of the corresponding extrapolation factors are in the range from 5% 
to 10 %. A double Gaussian fit is used to estimate the systematic errors of the extrapolation 
into the high y-region. The antiproton multiplicity normalized to (N w ), shown in the upper 
panel of Fig. 8, increases by nearly a factor 1.5 when going from mid-central Pb+Pb to 
inelastic p+p collisions and stays constant from mid-central to central Pb+Pb collisions. 
Although the properly weighted average of yields in p+p and n+n collisions would be the 
appropriate reference for this comparison, we consider the p+p midrapidity yield to be 
a good approximation. It is interesting to note that the antilambda multiplicity (from 
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FIG. 7: Rapidity density distributions per wounded nucleon are shown for protons at 40^4 GeV (top 
panel), protons at 158A GeV (middle panel), and antiprotons at 158^4 GeV (bottom panel) for five 
centralities in Pb+Pb collisions. The open symbols are obtained by reflection at mid-rapidity. The 
solid lines represent results of single Gaussian fits. For the statistical errors see Fig. 6. For the 
systematic errors see section IIIC. 

reference [21]), also shown in Fig. 8a, exhibits a centrality dependence similar to the one of 
the antiprotons. The lower panel of Fig. 8 shows the centrality dependence of the widths 
of the rapidity distributions for antiprotons and antilambdas (from reference [21]). No 
significant difference is observed between the values in different centrality bins of the Pb+Pb 
data. Only the width in p+p interactions seems to be slightly wider than those observed in 
Pb+Pb collisions. 

Net-proton rapidity distributions for all five centralities are obtained by subtracting the 
antiproton distributions from those of the protons. The result at 158A GeV for central 
collisions (bin CO) is compared to data published earlier [3] in Fig. 9. The differences 
between the two measurements can be traced back to different analysis methods. The earlier 
analysis used a method in which distributions of negatively charged particles were subtracted 
from those of the positively charged ones assuming the proton mass for all particles. The 



A 
V 



0.02 



CD 
DC 




100 200 300 400 



N 



w 



FIG. 8: The total multiplicities (per wounded nucleon) of p and A (upper panel) and the RMS y 
widths of the dn / dy distributions of p and A (lower panel) for five centralities as function of the 
number of wounded nucleons {N w ) in Pb+Pb collisions at 158A GeV. The A data are from [21]. 
Also shown are results on antiprotons obtained from p+p interactions [22]. Only statistical errors 
are shown (if larger than the symbol size). For the systematic errors see section IIIC. 

resulting distributions were corrected for the then unmeasured differences between 7r + and 
7r~, as well as K + and K~ yields by means of model calculations and detector simulations 



(see reference [3]). Here we identify the protons and antiprotons directly by means of their 
specific energy loss in the MTPCs. The new results should therefore be more reliable due to 
the smaller and better determined corrections leading to smaller systematic uncertainties. 
The two results are consistent within the systematic errors. 
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FIG. 9: The net-proton rapidity distribution at 158^4 GeV for centrality bin CO is compared with 
previously published data [3]. The full circles indicate results from the analysis presented in this 
paper, whereas the full squares show the previously published data. The open points result from 
reflection of the data points at midrapidity. The shaded bars represent the sytematic uncertainties. 
Error bars represent statistical errors. 



C. Model comparisons and conclusions 

Net-proton spectra are compared with HSD [23] and UrQMD-2.3 [24] model calculations. 
Before addressing possible differences between the model results and Pb+Pb collision data, 
we check how well the models reproduce the rapidity distribution of net-protons in p+p 
interactions. It is important to note here that the model calculations include also elastic 
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FIG. 10: Rapidity distribution of net-protons from p+p-collisions at 158 GeV/c beam momen- 
tum [22]. The solid line represents the HSD [23] and the dashed line the UrQMD [24] calculation. 

scattering whereas the experimental data represent inelastic interactions only. Fig. 10 reveals 
significant differences between the models and the experimental data. First, UrQMD [24] 
has more stopping than HSD [23] with a similar shape as the data, but a significantly 
higher yield at midrapidity. Second, HSD reproduces the midrapidity yield but fails to 
reproduce the shape. These shortcomings render comparison to Pb+Pb data questionable. 
The comparison is shown nevertheless in Fig. 11 for AO A GeV and in Fig. 12 for 158A GeV, 
in which the rapidity density distributions are divided by the average number of wounded 
nucleons (N w ). Note that antiprotons were not subtracted in the 40 A GeV data, since their 
yield is negligibly small compared to the proton yield (< 1% at midrapidity). Overall we 
observe fair agreement between the HSD model calculation and the data in the rapidity 
region covered by the NA49 measurements. Larger disagreements are seen for the UrQMD 
predictions. The differences can be attributed either to known differences between the 
models (see below) or to the deficiencies in describing proton distributions in elementary p+p 
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FIG. 11: The rapidity distributions scaled with 1/{N W ) of net-protons for five different central- 
izes at 40A GeV are shown together with results of HSD (top) [23] and UrQMD (bottom) [24] 
calculations. For the statistical errors see Fig. 6. For the systematic errors see section IIIC. 

interactions. The HSD [23] calculations describe the data always better than UrQMD [24] 
calculations. These significant differences in the rapidity distributions between HSD and 
UrQMD are due to a different definition of "formed" and "unformed" hadrons. In HSD a 
hadron is considered as "formed" only if the energy density (in the surrounding cell) drops 
below a critical value, which is taken to be 1 GeV/fm 3 in line with results from lattice QCD 
on the critical energy density for deconfinement. Otherwise the hadron is considered as 
"unformed" and unable to interact with other hadrons until the energy density drops to the 
critical value. This energy density criterion is not included in UrQMD which leads therefore 
to a substantial overestimate of the energy loss of participant nucleons in collisions in which 
high energy densities are reached [23]. 

In Pb+Pb collisions at 158 A GeV we observe a parabolic shape near to midrapidity (\y\ < 
1) which is similar to the one seen in inelastic p+p interactions as well as in HSD calculations. 
This can be explained by an exponential decrease of the proton yield towards midrapidity 
(according to exp[-(y — Ubeam)])-, which in turn results from the approximately flat probability 
density distribution as a function of Feynman x [1] . At 40A GeV the rapidity distributions 
are similar to those at 158A GeV for more peripheral collisions, while they develop a double 
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FIG. 12: The rapidity distributions scaled with l/(iV w ) of net-protons for five different centralities 
at 158,4 GeV are shown together with results of HSD (top) [23] and UrQMD (bottom) [24] 
calculations. For the statistical errors see Fig. 6. For the systematic errors see section IIIC. 

hump structure for central collisions. The normalized mid-rapidity yield increases by about 
25% at 40A GeV and 10% at 158A GeV from the most peripheral to the most central 
collisions in line with findings in earlier analyses of proton yields at mid-rapidity [5] and 
recent results on hyperons [21]. HSD calculations reproduce this trend quantitatively. 

The significant differences between the proton distributions at 40A GeV and 
158A GeV are at least partly due to the differences in experimental acceptances at the 
two energies. The change in shape which occurs away from midrapidity at 40A GeV is 
probably also present at the higher energy, however, in a range \y\ < (y max — 1.2) which 
is outside the NA49 acceptance. Finally, the integral of the normalized net proton distribu- 
tions over the measured region increases by 8% at A0A GeV and by 5% at 158A GeV with 
decreasing centrality. This is an indication that the spectrum of net protons contains nucle- 
ons not only from the nuclear overlap region specified in the Glauber model (which we call 
wounded nucleons) but also from the spectator remnants. Such an additional contribution 
to the proton spectra is confirmed in the results from the model calculations in which the 
coverage extends from target to beam rapidity (see Figs. 11 and 12). The origin of these 
"participating spectators" could be elastic and inelastic N+N spectator and meson+N spectator 



interactions (N stands for nucleon). Since there is no way to separate the "participating" 
spectator protons from the wounded nucleons in the experimental data, it is difficult to draw 
conclusions on nuclear stopping from the study of net proton spectral shapes in centrality 
selected Pb+Pb collisions. On the other hand this spectator contribution seems to be absent 
in the hyperon rapidity distributions [21], thus the e.g. meson+N speciator interactions are 
not violent enough to produce hyperons. 



V. SUMMARY 



The NA49 collaboration analysed proton and antiproton spectra in 40A GeV and 
158v4 GeV Pb+Pb reactions covering the 43.5% most central collisions. In the transverse 
mass spectra no strong variation with centrality is discernible, but the average tut in- 
creases by roughly 20% from peripheral to central collisions. The rapidity distributions at 
158^4 GeV have a common concave shape (Fig. 6b) which however gets shallower with in- 
creasing centrality, best seen if normalized to the number of participants as shown in Fig. 7. 
The rapidity distributions at 40A GeV on the contrary exhibit a strong centrality depen- 
dence starting with a concave or 'V shaped structure at large impact parameters, which 
turns into a symmetric double hump shape for more central collisions. The minimum at 
the center, which is characteristic of the concave shape, persists at all centralities. Sizeable 
contributions from non-participants in the rapidity range \y\ < (ybeam — 1-2) are observed, 
which are probably due to secondary interactions of produced particles in spectator matter. 
Thus rapidity loss analyses of net proton spectra in nucleus-nucleus collisions have to ac- 
count for this unwanted component of the spectrum. The midrapidity yield normalized to 
the number of participants varies only by 10% (5%) with centrality at 40A GeV (158A GeV) 
reaching its maximum for the most central collisions. As expected it decreases when going 
from 40 A GeV to 158 A GeV (by 25%). These finding may be compared to the correspond- 
ing data from AGS experiment E917 extracted from reference [2] using reference [25] to 
compute the number of participants. We find roughly a factor of two higher (normalized) 
yields at midrapidity with centrality variations of 3%, 5%, and 5% at 6A GeV, 8A GeV, and 
10. 8A GeV, respectively. Antiproton spectra could not be extracted at 40A GeV because 
of lack of statistics. At 158A GeV the normalized total (see Fig. 8) and midrapidity yields 
increase with impact parameter. Such a behavior is expected in case absorption plays a 



significant role or, equivalently, the baryon rich collision system approaches chemical equi- 
librium [20]. 

We chose to compare our findings to the transport model calculations from the HSD 
and UrQMD codes. Although both models reproduce rapidity density distributions in p+p 
collisions only with significant deviations as shown in Fig. 10, HSD gives a good description 
of the nuclear collision data at both energies (see upper rows in Figs. 11 and 12) in the 
region covered by NA49 measurements with deviations only at large rapidities. The UrQMD 
model calculations also fit the experimental distributions fairly well, although with larger 
discrepancies which are most pronounced at midrapidity and in central collisions. 

TABLE II: Proton rapidity densities dn/dy at five different centralities and 40 A GeV. The quoted 
errors are statistical. For the systematic errors see section IIIC. 





dn/dy(C0) 


dn/dy(Cl) 


dn/dy(C2) 


dn/dy(C3) 


dn/dy (C4) 


-0.02 < y < 0.18 


38.64±1.11 


30.69±0.79 


21.01±0.51 


12.69±0.33 


8.19±0.20 


0.18 < y < 0.38 


40.15±1.07 


30.63±0.75 


21.86±0.51 


13.26±0.34 


8.74±0.22 


0.38 < y < 0.58 


41.23±1.08 


32.93±0.70 


22.72±0.51 


14.49±0.34 


9.68±0.23 


0.58 < y < 0.78 


41.72±0.98 


34.24±0.74 


24.01±0.52 


15.48±0.35 


10.51±0.23 


0.78 < y < 0.98 


42.08±1.32 


33.65±0.72 


24.57±0.49 


16.64±0.37 


11.71±0.24 


0.98< y < 1.18 


40.11±1.19 


32.91±0.85 


24.58±0.50 


17.48±0.38 


12.71±0.28 


1.18 < y < 1.38 


36.64±0.98 


31.29±0.73 


24.84±0.53 


18.70±0.42 


14.06±0.33 


1.38 < y < 1.58 


31.63±1.09 


29.16±0.76 


24.94±0.56 


20.69±0.51 


15.69±0.48 



TABLE III: Proton rapidity densities dn/dy at five different centralities and 158^4 GeV. The 
quoted errors are statistical. For the systematic errors see section IIIC. 





dn/dy(C0) 


dn/dy(Cl) 


dn/dy(C2) 


dn/dy(C3) 


dn/dy(C4) 


-0.52 < 


y 


< 


-0.32 


25.63±1.41 


22 


17±0 


74 


14.36±0.43 


9.28±0.30 


5.84±0.21 


-0.32 < 


y 


< 


-0.12 


27.49±1.24 


22 


32±0 


61 


15.14±0.38 


9.46±0.25 


6.03±0.18 


-0.12 < 


y 


< 0.08 


27.51±1.09 


20 


97±0 


52 


14.79±0.33 


9.46±0.23 


6.07±0.19 


0.08 < 


y 


< 


0.28 


28.07±1.01 


21 


81±0 


51 


14.96±0.31 


9.34±0.21 


5.98±0.16 


0.28 < 


y 


< 


0.48 


27.93±0.94 


21 


74±0 


52 


15.35±0.30 


9.64±0.21 


5.96±0.15 


0.48 < 


y 


< 


0.68 


29.20±0.89 


22 


91±0 


50 


16.27±0.30 


9.92±0.23 


6.39±0.16 


0.68 < 


y 


< 


0.88 


31.30±0.88 


23 


96±0 


49 


17.4±0.31 


10.97±0.23 


7.05±0.17 


0.88 < 


y 


< 


1.08 


30.51±0.85 


25 


10±0 


49 


17.87±0.32 


11.53±0.24 


7.55±0.19 


1.08 < 


y 


< 


1.28 


32.35±1.05 


26 


88±0 


57 


19.22±0.39 


12.25±0.26 


8.43±0.20 


1.28 < 


y 


< 


1.48 


38.84±3.43 


28 


82±0 


58 


21.46±0.43 


13.63±0.28 


9.52±0.22 


1.48 < 


y 





1.68 


339.35±2.00 


30 


30±1 


55 


20.13±2.30 


16.38±0.98 


10.73±1.31 



TABLE IV: Antiproton rapidity densities dn/dy at five different centralities and 158^4 GeV. The 
quoted errors are statistical. For the systematic errors see section IIIC. 





dn/dy(C0) 


dn/dy(Cl) 


dn/dy(C2) 


dn/dy(C3) 


dn/dy(C4) 


-0.52 < y < -0.12 


1.61±0.13 


1.37±0.09 


0.92±0.06 


0.78±0.05 


0.51±0.04 


-0.12 < y < 0.28 


1.75±0.12 


1.37±0.08 


1.06±0.05 


0.80±0.04 


0.55±0.03 


0.28 < y < 0.68 


1.76±0.11 


1.43±0.07 


0.96±0.05 


0.70±0.04 


0.49±0.03 


0.68 < y < 1.08 


1.31±0.10 


1.11±0.07 


0.80±0.05 


0.54±0.038 


0.36±0.03 


1.08 < y < 1.48 


0.60±0.10 


0.57±0.09 


0.37±0.05 


0.21±0.03 


0.31±0.06 
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